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Types  of  distributed  radar  surveillance  systems  range  from  synthetic 
aperture  radars  (SAR's)1  to  netted  systems  of  monostatic  radars.2  Two 
obvious  discriminants  for  classifying  this  range  of  types  are  the  amount 
of  system  coherence  in  the  sensor  signal  processing  and  the  degree  of 
autonomy  allowed  the  various  elements.  The  SAR,  which  is  certainly  a 
distributed  system  (as  a  function  of  time ) ,  processes  all  sensor  signals 
in  a  single  coherent  processor;  however,  the  typical  SAR  can  be  viewed 
as  a  rigidly  constrained  (time)  sequence  of  locations  of  monostatic 
radars.  Thus,  the  "elements”  of  a  SAR  possess  zero  autonomy.  The 
netted  system  of  monostatic  radars,  wherein  each  radar  ignores  echo 
signals  generated  by  all  other  radars,  has  no  system  coherence;  in  fact, 
the  system  discards  or  ignores  potentially  useful  signals  but  allows  an 
extreme  autonomy  of  its  various  elements.  An  interesting  feature  of 
these  two  systems  is  the  fact  that  the  map  produced  by  a  SAR  is  self- 
registered  and  relatively  independent  of  the  path  followed  by  the  SAR; 
thus,  only  the  relative  precision  of  successive  SAR  locations  is 
constrained.  On  the  other  hand,  the  netted  radar  system  requires 
knowledge  with  absolute  accuracy  of  the  locations  of  all  its  elements. 

Between  these  two  extremes  lie  a  number  of  interesting  distributed 
surveillance  systems  whose  configurations  are  pertinent  to  various 
tactical  and  strategic  scenarios.  Especially  interesting  are  systems 
that  can  obtain  precise  target  coordinates  without  the  use  of  very 
narrow  antenna  beams,  thus  allowing  some  reduction  in  antenna  size. 
Some  of  these  systems  are  described  and  analyzed,  in  which  time-interval 
measurements  alone  are  sufficient  to  allow  complete  resolution  of  their 
geometry.  While  we  intend  only  to  address  geometry,  nongeometric  radar 
facts  will  be  brought  to  bear  as  needed  to  keep  the  systems  realistic. 
For  the  systems  studied  here,  it  will  be  seen  that  precise,  inaccurate 
(nonsynchronized)  clocks  suffice;  a  single  benchmark  (accurate  location 
of  an  element  in  map  coordinates)  is  adequate;  and  system  coherence  can 
be  traded  for  data  redundance,  just  as  coherent/incoherent  signal 
processing  can  be  traded  off  in  a  monostatic  radar.  An  optimal 
exploitation  of  data  redundance  is  developed  for  each  system. 


1John  J .  Kovaly,  Synthetic  Aperture  Radar ,  Aztech  House  (1976). 

^Netted  Radar  Program,  Volume  lit  TPS-5X  Ground  Surveillance  Radar, 
MIT  Lincoln  Laboratory,  contract  F19628-80-C-0002  (30  September  1981). 


2.  NOTATIONS  AND  DEFINITIONS  OF  PARAMETERS 

The  following  symbols  are  used  in  this  report. 

(1)  x,  y,  z  Conventional  Cartesian  coordinates 

(2)  j,  k,  a  Integer  subscripts  identifying,  respectively,  (1)  a 

distinct  source  at  positions  Pj  ■  Xj,  yu,  z • ;  (2)  a 
distinct  sensor  (receiver)  at  positions  =  xk, 
yk,  zk,  and  (3)  a  distinct  target  at  positions  Pa  = 
xa,  ya,  za.  Integers  m,  n,  and  q  are  defined  such 
that  0  <  j  <m,  0  <  k  <  n,  and  0  <  a  <  q. 

(3)  C  A  central  processing  station  at  position  Pc  = 

(xc,yc,zc)  =  (0,0,0).  A  (one-way)  communications 
link  exists  (where  appropriate)  from  every  sensor 
or  transmitter  position  Pj  to  the  central  proc¬ 
essing  station. 

A  time  interval  multiplied  by  the  speed  of  light 
(in  meters  per  second)  and  read  out  in  the  system 
as  slant  range  from  position  Pu  to  position  Py. 

A  vector  of  magnitude  |Rft  |  *  RUv  *  Tuv  from 
position  Pu  to  position  Py. 

(6)  A  Absolute  time  measured  from  the  beginning  of  the 

universe.  Thus,  the  counter  of  a  conventional 
cyclic  clock  will  indicate  A(modulo  M)  +  B,  where  B 
is  a  bias  in  the  clock  count  that  has  developed  due 
to  long-term  drift,  and  M  equals  a  full-scale  (plus 
one  count)  state  of  the  counter.  Parameters  A,  M, 
and  B  are  measured  in  meters  and  represent  the 
distance  a  light  pulse  would  travel  during  the 
clock  periods  represented  by  A,  M,  and  B. 

(7)  N  A  number  indicated  by  a  clock,  i.e.,  the  state  of 

the  clock -counter.  A  clock  exists  at  each  element 
of  the  system. 


3.  LOCATING  SOURCES 

In  this  and  succeeding  analyses  of  distributed  surveillance  systems, 
functional  relations  between  various  vectors,  coordinates,  and  other 
parameters  are  presented  in  equation  form.  These  equations  are  used  to 


(4)  Tuy 

(5)  *uv 


t 


(1) 


where  R.^  is  the  distance  in  meters  separating  source  j  from  sensor  k, 
and  RjcJi8  the  distance  separating  source  j  from  C. 

Equation  (1)  will  be  used  to  derive  a  set  of  n  linear  equations  in 
unknowns  ,  y  ■ ,  z  • ,  and  RjC«  In  order  to  accomplish  this,  a  vector  V 
with  elements  is  defined  such  that 

vi  -  -  fjc)2  *  D1 

V2  ■  (T:a  -  Tjc)2  -  °2 

•  •  •  (2) 


Vk  -  (Tjk  -  Tjc)2  -  4 

Vn  «  (Tjn  -  Tjc)2  -  oj  . 


The  differences  of  adjacent  vector  elements  are 

V,  -  V2  -  (TJt  -  Tje)2  -  (TJ2  -  TJc)2  , 

V2  *  V3  *  (Tj2  *  Tjc)2  '  (TJ3  -  Tj=)2  ’ 

•  •  • 

.  .  .  (3) 

•  •  • 

Vk  “  Vk+1  "  (Tjk  “  Tjc)2  "  ^Tj  (k+1 )  "  Tjc^2 

•  •  • 

•  •  • 

•  •  • 

v  ,  -  V  -  (T.  .  ..  -  T.  )2  -  (T.  -  T.  )2 

n-1  n  i  j(n-1)  jc'  ^  jn  jo' 
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*• 


0  0-  '<0.-  0  -0  '  -'•0 


Expanding  the  k"1  element  of  (3)  gives 


^  “  'Vk+I  -  (^k  -  (k+1 ))  “  2Tjc(Tjk  "  Tj  (k+1 ))  *  (4) 


which,  in  accordance  with  definition  5,  becomes 

2  2 

Vk  -  Vfc+1  =  Rjk  -  Rj(k+1)  -  2Rjc(Tjk  -  Tj(k+1))  •  (5) 

It  follows  from  the  geometry  in  figure  1  that 

*jk  =  (xj  “  xjc)2  +  (yj  -  yk)2  +  (zj  -  zk)2  (6) 

and 

Rjtt+i)  “  (Xj  -  xk+1)2  +  (yj  -  yjc+i)2  +  (Zj  -  Zk+1)2  .  (7) 

Also,  since  xc  =  yc  *  zc  »  0, 

2  2  2  2 

^k  *  xk  +  yk  +  zk  (8) 

and 

Kcfk+I)  *  *k+1  +  yJ+1  +  zk+1  •  O) 

Furthermore ,  by  equation  ( 1 ) , 


T  -  T  ■  D  -  D 

jk  j (k+1 )  k  k+1 


(10) 


Combining  equations  (5)  to  (10)  and  collecting  terms  in  the  unknown 
quantities  Xj,  yj,  Zj,  and  RjC  gives 


2(*k+i  -  xk)*j +  2('rk+i  -  »khj  *  2(*k+i  -  zk)2j 


(ID 


+  2(0)5+!  -  Dj5)Rjc  -  [Vjt  -  Vk+T  +  (R|(k+D  -  R^k)] 


'iV_-.^: *  J.-.J.’’  -r*- i  -  V>t>.li  **, 


The  known  quantities  in  equation  (11)  are  the  expressions  on  the 
right  of  the  equality  sign  and  the  coefficients  of  Xj,  y j ,  z- ,  and  Rjc 
on  the  left.  The  quantities  and  in  the  Rjc  coefficient  are 

determined  by  direct  measurement.  J 

Applying  these  results  to  all  elements  of  equation  (3)  gives  the 
following  equations: 

2(x2  -  x1)xj  +  2(y2  -  yjy.  +  2(z2  -  Zl)Zj  +  2(D2  -  D^R^ 


Vi  -  V2  +  R22  -  Rci 


-  *  2(yk+i  -  yJyj  +  2(Zk+1  "  +  2(Dk 


2  2 

vk  "  Vj,+i  +  Rj(Jt+1)  “  Rck  f 


m 


2K  -  V.l'j  +  2fyn  -  Vl)yJ  *  2(*»  -  V,)2j  *  2(D„  -  D„-1>Rj 


vn-1  “  vn  +  Rcn  “  Hc(n-1) 


Equations  (12)  constitute  a  set  of  n  -  1  simultaneous  linear 
equations  in  the  four  unknowns  x . ,  y . ,  z j ,  and  Rjc*  The  equations  are 
solvable  if  n  >  5.  Thus,  the  position  Pi  of  the  jth  source  can  now  be 
found  each  time  it  emits  a  pulse.  Therefore  Pj  can  vary  from  pulse  to 
pulse,  either  by  variation  of  the  position  of  a  single  source,  or  by 
variation  of  the  index,  j,  of  a  number  of  sources. 
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2 


(x2  -  xl) 

(*2  -  yi) 

(z2  -  zl) 

(°2  -  Dl) 

Xj 

• 

• 

• 

• 

yj 

(xk+1  -  xk) 

+ 

i 

(zk+1  -  zk) 

(°k+1  *  °k) 

zj 

• 

(*n  -  *n-l) 

r~ 

1 

• 

1 

• 

(zn  “  zn-l) 

• 

(Dn  -  Dn-l) 

Rjc 

( v,  -  V2  +  r|2  -  *4l) 

•  •  • 

2  2 

(Vk  -  Vk+1  +  RcOt+i)  -  Rck) 
•  •  • 

(vn-1  “  vn  +  Rcn  “  ^(n-l)) 


(13) 


If  n  exceeds  4  in  equation  (13),  an  overconstrained  set  of  equations 
exists,  and  it  is  almost  impossible  to  conceive  of  a  single  vector  (Xj, 
y, ,  2, ,  r.c)  that  when  multiplied  by  the  row  vectors  of  the  matrix  in 
(13)  wouldgive  precisely  the  measured  values  in  the  vector  on  the  right 
side  of  (13).  Clearly,  some  estimated  vector  (x^,  y^,  z^,  Rjc)  is 
required  that  fits  the  v'  le  data  set  in  the  measurement  vector. 
Equation  (13)  can  be  written  in  the  compact  symbolic  matrix  format 


2[A](J]  =  [B] 


or 


2AJ  =  B 


The  solution  requires  a  guess  vector,  G,  such  that 


2AG  =  C 


However,  this  produces  an  error  vector 


E  =  C  -  b  , 
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the  elements  of  which  must  be  minimized  in  some  fashion.  Gauss  solved 
this  problem*  using  the  criterion  that  the  sum  of  squared  elements  of  E 
is  minimized.  To  see  how  this  was  done*  note  that 


E2  =  E_E  , 
T 


where  E^  is  the  transpose  of  the  vector  E.  Thus, 


ETE  =  (2AG  -  B )T( 2AG  -  B) 


Differentiating  partially  with  respect  to  G  gives 

*(V0 

— —  =  2A_ ( 2AG  -  B)  +  (2AG  -  B )  2A 
3G  T  T 


=  2At(2AG  -  B)  +  2At(2AG  -  B) 


=  8A^pAG  -  4AtB  . 


Equating  this  derivative  to  zero  provides  the  desired  result: 


2AtAG  *  A<j3  . 


The  matrix  A^A  is  a  square*  4x4  nonsingular  matrix.  Thus,  equation 
(14)  provides  a  unique  solution  for  G,  independent  of  the  number  of 
sensors.  The  solution  for  G  may  be  inserted  into  equation  (13),  and  the 
guess  vector  2AG  =  C  allows  the  evaluation  of  C  -  B  as  a  check  on  the 
overall  system.  If  it  is  found  that  two  adjacent  elements  of  C  -  B  are 
unduly  large*  then  it  may  be  concluded  that  the  sensor  with  clock 
readings  common  to  these  two  elements  has  been  mislocated*  or  has  a 
clock  malfunction*  etc. 

Further  discussion  of  the  matrix  [A]  is  appropriate  at  this  point. 
Hie  first  three  column  vectors  of  [A]  represent  static  data  so  long  as 
no  sensors  are  moved;  thus,  only  the  fourth  column  vector  of  [A] 
contains  real-time  ( pulse -to-pulse )  information.  If  the  distance 
between  each  sensor  and  C  could  somehow  be  directly  measured,  then  [A] 
would  become  a  three-column  matrix,  and  A^A,  along  with  its  inverse, 
could  be  precomputed  just  once  for  all  sources  or  source  positions,  thus 


'■■-ju-  J~£  _  .  M 


reducing  the  real-time  computation  load  by  a  considerable  factor.  Hie 
implications  of  this  fact  will  be  developed  at  a  later  point  in  the 
discussion. 


4.  LOCATING  A  SINGLE  TARGET 

Before  we  discuss  system-target  geometry,  we  should  point  out  that, 
in  general,  we  are  no  longer  discussing  single-pulse  timing.  If  there 
were  no  clutter,  and  if  each  target  were  large  enough  in  radar  cross 
section  to  bring  its  echo  clear  of  thermal  noise  in  the  receivers,  then 
we  could  proceed  pulse  by  pulse.  Since  clutter  is  present  and  targets 
are  small,  it  is  required  that  each  sensor  have  some  target-detection 
mechanism  that  separates  meaningful  targets  from  clutter,  reporting  the 
average  clock  reading  that  corresponds  to  each  target's  position,  Pft. 
Clearly,  if  there  is  to  be  motion  of  the  source (s),  the  requirement  is 
that  the  target-detection  interval  at  each  sensor  be  short  enough  to 
make  source  motion  during  each  interval  negligible. 

Assume  now  that  (a)  the  system  is  as  described  in  the  previous 
section  and  (b)  the  location  of  the  source(s)  is  now  known  (see  fig. 
2).  A  single  target.  A,  is  added  at  position  Pa  (xa»Ya»za) 
unobstructed  line  of  sight:  from  the  target  to  the  source  (s)  and 
sensor(s).  Assume  further  that  each  sensor  detects  the  target  after 
receiving  a  sequence  of  echos  due  to  one  sequence  of  illumination  pulses 
from  source  j .  Each  sensor  then  sends  a  clock  reading  to  the  central 
processor,  C,  that  represents  the  total  distance  from  Pj  to  P^  via  Pa. 


This  distance  is  Rja  +  Rftk.  At  position  C,  each  reading  is  squared  and 
a  vector,  B^,  of  successive  differences  of  these  squares  is  formed. 
Thus 


\  *  K.  ♦  R, 


ja  a(k+1) 


)2  -  (RJ,  *  Rak)2 


=  2Rja(Ra(k+1)  “  Rak)  +  Ra(k+1)  ”  Rak  • 


The  term  2Rja(Ra(k+1)  "  Rak)  comPri-sea  the  unknown  R.  and  twice  the 
difference  Between  two  knowns,  Ra (K-l-1 )  and  Rak*  difference  in 

terms  of  attendant  coordinates  at  positions  Pa,  Pk,  and  Pk4,  is 
obtainable  from 


2.  2 

Ra(k+1 )  -  Rak  =  (xa  -  xk+l)2  +  faa  -  Vk+l)2  +  (2a  *  2k+l)2 

\)2  ♦  (ya  -  yk)2  +  (■.  -  *k)2] 

In  this  equation,  we  may  substitute 


2  2  f  2  2  2  2  2  2 

Rc (k+1 )  -  Rck  =  (xk+1  +  yk+1  +  2k+l)  “  (xk  +  Yk  +  2k)  » 


^  2  2  jl 

Ra(k+1)  -  Rak  -  Rc(k+1)  “  Rck  “  2(xk+1  ~  xk)xa 


-  2(yk».  -  yk)y.  -  2Kt1  -  \K  ■ 

Combining  equations  (15)  and  (16)  and  collecting  terms  in  all  unknowns 
xa'  ya'  2a'  and  Rja  9ives 

2K*1  -  "k^a  +  2(*k.,  -  +  2(v,  -  2kK  -  2^Ra(k+l )  -  Bak^Rja 

1  2  (,7, 

*  Rc  ( k+1  )  “  Rck  “  Bk  • 


ti si.y; r  -mt, -m* 


Equation  (17)  can  be  written  successively  for  the  cases  K  >  1,  2,  ...k 
...  n,  thus  yielding  the  matrix  equation 


(x2  -  *,)  (y2  -  yj  (*2  ■  zl)  (Ra1  ■  Ra2) 


(xk+1  -  xk)  (yk+1  -  yk)  (zk+1  -  zk)  (Rak  -  Ra(k+1)) 


(^Si  -  xn-l)  (yn  ■  yn-l)  (zn  “  zn-l)  (Ra(n-1)  "  Ran) 


Rja 


(rc2  -  *4l  -  ®l) 


e  e 


CRc(k+1 )  “  Rck  *  Bk) 


(18) 


I  (®cn  ~  ®c(n-1)  "  ®n-l)  | 

This  matrix  equation  can  be  solved  by  the  methods  used  previously  in 
solving  equation  (13). 


Note  here  that  the  matrix  of  known  values  in  equation  (18)  differs 
from  the  matrix  of  (13)  only  in  the  fourth  column.  Also,  if  the  range 
from  the  source  to  the  target,  R . a,  were  known  by  direct  measurement, 
then,  analogous  to  the  source-location  task,  the  real-time  computation 
load  would  diminish  considerably. 


5.  SINGLE  AUTONOMOUS  SENSOR 


Harking  back  to  the  system  described  in  section  4,  a  single  roving 
sensor  is  now  postulated  that  contains  all  the  computing  power  of r 


central  processor  C  and  also  "knows”  its  own  location  in  absolute  map 
coordinates  at  all  times  (see  fig.  3).  If  the  successive  (known) 
positions  of  this  sensor  are  labeled  x^,  yk,  (1  <  k  <  n),  then  after 
all  n  positions  have  been  occupied  by  this  roving  sensor,  a  computation 
can  be  made,  as  in  section  3,  to  determine  the  relative  locations  of  the 
sources.  In  this  case  T. ^  «  Df  is  substituted  for  -  T.  ■  D^. 
Since  Tjc  =  Rjc  no  longer  enters  into  the  computation,  tne  fourth  column 
of  the  matrix3 in  equation  (13)  is  not  required  in  the  computations. 
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Figure  3.  Distributed  surveillance  system 
(single  autonomous  sensor  ) • 


Having  thus  determined  the  matrix  of  relative  source  locations,  we 
may  assume  some  favorable  (known)  vantage  position  for  seeking  out 
targets.  Once  again,  if  a  single  target  is  assumed,  then  for  each 
source  j  a  separate  clock  reading  of  Tja  +  TaJt  ■  Rja  +  ^  ■  Fj  will  be 
obtained.  Computing  much  as  before,  we  find 


Fj+1  -  “  R*j+1)a  "  Rja  +  2Rak(R(j+1)a  “  Rja) 


R(j+1)a  ”  Rja  “  xj+1  "  xj  “  2xa(xj+1  "  xj) 


2  2 

+  yj+i  -  yj  -  2ya(yj+i  -  yj) 
2  2 

+  Zj+1  ”  zj  *  2za(zj+1  “  zj) 


Here,  the  matrix  equation  requiring  solution  is 

(x2  -  xt)  (y2  -  Yi)  (z2  -  zl)  (R1a  "  R2a) 

e  e  e  e 

e  e  e  e 

2  (xj+i  -  xj)  (yj+i  -  (*j+i  -  zj)  (*j«  -  R( j+i  >a) 

e  e  e  e 

e  e  •  • 

(*n  -  *n-l)  (yn  -  yn-l)  (zn  -  xn-l)  (R(n-1 )a  "  «na) 

(p!  -  P?  ♦  R 1 ( 2 ) ) 

e  e  e 

•  e  • 

(y§+i  -  ♦  Rj<j+i)) 

•  •  e 

•  #  e 

(*n  "  *n-1  +  R(n-1)n) 


The  solution  is  carried  out  as  in  the  previous  systems.  Note  that  each 
of  n  independent  autonomous  sensors  can  independently  exploit  a  field  of 
sources  in  this  manner,  and  the  additional  cost  of  replicating  the 
computing  power  of  C — and  its  ability  to  determine  its  own  location — may 
well  be  overcome  by  the  fact  that  these  sensors  can  be  utterly  passive. 


6.  REMOVAL  OF  CLOCK  BIAS 

The  initial  deployment  of  a  distributed  system  might  be  accomplished 
in  a  number  of  ways.  Conceptually,  the  easiest  deployment  would  use 
benchmarks  that  have  been  accurately  surveyed  on  the  ground;  however,  a 
system  that  required  such  a  survey  would  be  severely  limited  in  its 
range  of  applications.  There  is  a  technique  that  allows  both  self- 
established  deployment  and  continual  monitoring  of  incremental  modifi¬ 
cations  to  accommodate  the  expected  flux  of  a  battlefield  situation. 
Since  this  technique  can  establish  only  the  relative  locations  of  the 
system  elements,  one  is  required,  as  might  be  expected,  to  know  the 
absolute  location  (map  coordinates)  and  orientation  (north)  of  at  least 
one  element  of  the  system.  For  ease  of  discussion,  assume  that  the 
central  processing  element  (Pc)  is  the  one  whose  map  coordinates  and 
orientation  are  known.  A  high-resolution  three-dimensional  radar  is 
placed  at  C,  and  used  to  determine  the  pulse  source  positions,  P j . 

Consider  m  pulse  source  positions  Pj,  0  <  j  <  m;  receiver  positions 
0  <  k  <  n;  and  a  central  processing  position"  Pc«”  Clocks  exist  at  all 
positions  with  the  following  characteristics: 

(a)  All  clocks  have  "identical"  crystal -control led  frequency 
sources . 

(b)  Every  clock  has  a  b-bit  counter,  automatically  recycling  to  a 
zero  count  after  each  2 (exp  b)  cycles  of  its  frequency  source;  thus,  the 
clock  counters  have  modulus  M  =  2 (exp  b) 

(c)  Each  clock  has  a  bias  B: 

clock  at  Pj  has  bias  B j , 
clock  at  P^  has  bias  B^, 
clock  at  Pc  has  bias  Bc. 

Assume  that  all  sources  have  the  same  interpulse  interval,  and  that 
this  interval  is  equal  to  the  modulus  M  of  the  clock  counters*  Modulus 
M  is  selected  large  enough  to  avoid  ambiguous  ranging  and  troublesome 
second -time -around  echos.  Finally,  assume  that  each  receiver  can 
uniquely  determine  the  identity  of  the  source  of  each  pulse  it  receives. 


•1 

'i 


will  be  noted  as  Tj 
C  and  that  Tjc  "*j 


a8  ♦  ®k  "  B1  "  "ki* 

,  -  Rj_,  we  have,  for  ea< 


Assuming  that  all  Pj  are  known  at 
ch  pulse  at  C, 


*cj  "  Rjc  "  Bc  "  Bj 


It  is  assuaed  that  the  numbers  as  measured  at  the  receivers  are 
somehow  transmitted  to  C  (e.g.,  by  laser/  line  of  sight,  telephone 
lines,  and  so  on),  not  necessarily  in  real  time. 

This  allows  the  further  computation 

Nkj  -  ("Cj  -  Rjc)  “  Tjk  +  ®k  '  Bc  “  ^k  *  Bc  +  ®k  •  (21) 

An  examination  of  figure  4  shows  that,  so  long  as  no  single  time 
interval  exceeds  M,  and  so  long  as  all  time  computations  are  performed 
modulo  N,  these  timing  equations  are  correct,  regardless  of  the 
magnitudes  of  the  biases.  At  C,  the  analysis  proceeds  with  the 
determination  of  the  relative  locations  of  all  receivers  and  the  clock 
biases.  The  computations  are 


m 


L  ■  (Rjk  +  ®k  **  ®c)2  "  Rjk  +  B*  +  B2  +  2B)cRj]c  -  2BcB]c  -  2BcRj]c  , 


M  *  (R(j+1)k  +  “  Bc)2  *  R?j+1)k  +  ®k  +  Bc  +  2V*(j+1)k  "  2BcBk 


.2  . 


'  2BcR(j+1)k  ' 


N  -  (Rj (k+1 )  +  ®k+1  -  Bc)2  *  Rj (k+1 )  +  Bk+I  +  Bc 
+  2Bk+1Rj(k+1)  “  2BcBk+1  “  2BcRj(k+1)  , 


2  2  2 
p  ■  (R( j+1 )  (k+1 )  +  ®k+1  “  Bc)2  "  R(j+1)(k+1)  +  ®k+1  ♦  Bc 

+  2Bk+1R( j+1 ) (k+1 )  "  2BcBk+1  ”  2BcR( j+1 ) (k+1 )  ' 


M  -  L 


R(j+1)k  ”  Rjk  +  2(®k  -  BcJ(R(j+1)k  “  Rjk) 


, 


(26) 


and 

2  2 
P  -  L  -  R(j+i)(k+l)  -  RjCk+l) 

+  2(®k+1  "  BJ(R(j+1)(k+l)  "  Rj(k+1)) 


(27) 


Heref  the  coefficients  of  \  -  Bc  and  Bjt+^  -  Bc  are  known  via  direct 
measurement,  for 


R(j+l)k  "  Bc  +  ®k  _  (Rjk  “  Bc  +  ®k)  “  R(j+1)k  “  Rjk 

and 

R( j+1 ) (k+1 )  “  Bc  +  Bk+1  ”  ^Rj (k+1 )  “  Bc  +  ®k+1^ 

"  R(j+1)(k+1>  “  Rj(k+1)  * 


Moreover , 

[P  -  N]  -  [M  -  Ll  - 


2(xk  “  *k+1^xj+1  ”  xj)  +  2(yk  ’  yk-t1^yj+1  “  yj) 

+  2(zk  -  «k+1)(*j+1  -  *j)  +  2(®ic+i  ”  Bc^R(j+1)(k+l)  “  Rj  (k+1  )^ 
+  2(Bc  *  Bk)(R(j+1,k  -  Rjk)  *  ( 


The  left  side  of  this  equation  comprises  terns  that  are  computed 
solely  from  the  various  tine  measurements  reported  by  the  system.  The 
right  side  contains  two  clock  bias  terms,  each  with  a  known  coefficient, 
and  the  products  of  like  Cartesian  coordinates  of  source-to-eource 
distances  with  sensor -to -sensor  distances.  Clearly,  if  we  know  five  or 
more  of  either  set  of  distances,  we  have  a  solvable  set  of  equations  for 
the  unknown  distance  intervals. 
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The  foregoing  technique  can  be  readily  applied  to  the  case  of  a 
single  roving  sensor  that  contains  all  the  computational  power  of  the 
system*  if  we  require  that  this  sensor  have  knowledge  of  its  own 
position  in  map  coordinates  at  all  times*  and  constrain  the  sources  to 
be  immobile  at  (initially)  unknown  positions  Pj. 

This  sensor  will  record  its  own  clock  states  Nj^  at  a  sequence  of 
known  positions  P^  (xjc*y]c*z]c).  As  before*  each  source  emits  a  pulse  at 


Nj  ■  0  and  at  A 
and  the  state  wi 


Vj  ■  gM  —  Bj.  These  are  received  at  A^j 
ill  be  recorded  as 


gM  -  Bj  +  Tjk, 


“  Tjk  -  Bj  +  ®k 


Recall  now  that  is  the  bias  of  the  clock  in  the  single  roving  sensor; 
thus,  it  will  not  vary  as  location  Pk  is  varied.  Substituting  Bc  for  B^ 
in  this  equation  gives 


\j  “  Tjk  +  Bj  +  Bc 


R -  B  +  B 
jk  j  c 


This  equation  is  of  the  same  form  as  equation  (21)  except  for  the  change 
in  sign  of  Bc*  the  fact  that  -Bj  is  involved  rather  than  iB^,  and  the 
fact  that  the  P^  rather  than  the  Pj  are  now  known.  It  should  be  clear 
that  the  analysis  can  proceed  as  with  equation  (19)*  with  known  and 
unknown  position  coordinates  interchanged*  and  with  biases  Bj  rather 
than  B^. 

Having  determined  the  relative  positions  of  the  sources*  and  the 
biases  of  their  clocks  relative  to  its  own  clock*  we  may  program  the 
single  autonomous  sensor  to  continually  monitor  the  direct-path  signals 
from  these  sources  to  correct  for  further  changes  in  any  clock  biases. 
Similarly,  clock  biases  in  the  aultisensor  arrangement  with  a  central 
clock  can  be  continually  monitored  and  corrected. 


7.  REPORT  SORTING  FOR  MULTIPLE-OBJECT  DETECTIONS 


In  the  systems  so  far  discussed*  angular  measurements  have  been 
unnecessary.  Only  the  best  and  most  easily  attainable  aspect  of 
electromagnetic  surveillance  techniques  have  been  used  to  achieve  the 
accuracy  of  position  location  desired.  There  is  at  least  one 
application  that  requires  only  this  measurement  of  time;  range 


*- 
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instrumentation  for  tracking  the  location  of  one  or  more  cooperative 
sources.  The  multiple  sensor  case  is  ideal  for  this  task,  if  at  least 
one  nonmoving  source  is  retained  to  allow  tracking  of  the  sensors'  clock 
biases . 

The  more  general  application,  where  noncooperative  moving  objects 
must  be  detected  and  located,  requires  angular  measurements  to  some 
degree  for  two  reasons:  the  first  is  a  practical  radar  consideration 
involving  clutter  levels,  and  the  second  is  the  report-sorting 
problem.  Figure  5  depicts  the  dilemma  of  any  of  our  systems  when 
confronted  by  more  than  one  detection.  Assume  objects  A  and  B  lie  on 
isochrons  for  sensors  k  and  k  +  1  as  shown.  Each  sensor  will  report  two 
detections  at  clock  states  corresponding  to  the  appropriate  isochron. 
The  central  processor  can  now  associate  the  reported  clock  states  in  a 
fashion  that  indicates  objects  at  C  and  D,  as  well  as  at  A  and  B. 
Additional  sensors  will,  of  course,  provide  additional  data  in  more  than 
sufficient  quantity  to  resolve  the  simple  ambiguity  indicated  in  figure 
5,  but  consider  the  processing  load.  If  m  sensors  have  each  reported  n 
distinct  detections,  the  central  processor  has  nm  possible  combinations 
to  sort  through. 


ISOCHRONS  FOR  k  +  1 


ISOCHRONS  FOR  k 


Figure  5.  Multiple  target  detection  dilemma 


Figure  5  shows  sensors  k  and  k  +  1  with  multiple  simultaneous 
antenna  beams.  In  the  situation  depicted,  these  antenna  beams  resolve 
the  ambiguity,  for  each  sensor  now  reports  both  the  clock  state  and 
receiving  beam  identity  for  each  detection.  We  note  in  passing  that  the 
antenna  beam  orientation  of  each  sensor  is  approximately  deducible  at 
the  central  processor,  since  the  direct-path  geometry  is  known.  Since 
angle  measurement  is  not  being  used  to  establish  positions  of  objects, 
but  only  to  reduce  ambiguities,  an  approximate  knowledge  of  sensor 
orientation  is  adequate. 

The  architecture  of  the  report-sorting  mechanism  is  strongly 
affected  by  certain  practical  system  considerations.  Each  sensor  will 
have  as  few  antenna  beams  as  possible,  so  that  the  computational 
complexity  of  the  antenna  beam  former  is  minimized.  In  most  cases,  for 
objects  on  the  ground,  clutter-level  considerations  will  dictate  the 
irreducible  minimum  of  antenna  beams  required.  One  parameter  of  the 
system  that  will  be  fixed  by  the  choice  of  angular  resolution  is  the 
total  integration  time  for  detection  at  each  sensor.  The  maximum 
autonomy  of  sensors  is  desired;  therefore,  it  will  be  assumed  that  the 
integration  epochs  at  the  sensors  are  not  mutually  synchronized.  This 
means  that  the  central  processor  can  never  be  assured  that  it  has  a 
complete  static  data  set  for  any  given  object;  moreover,  the  elements  of 
the  data  set  for  a  given  object  will  have  an  age  variation  on  the  order 
of  the  integration  interval,  and,  of  course,  many  objects  will  not  be 
detected  by  all  sensors.  Finally,  for  full  exploitation  of  the  system 
possibilities,  each  sensor  must  have  a  reasonably  high  false-detection 
probability.  To  deal  with  these  and  other  complicating  factors  to  be 
elucidated,  the  surveillance  area  shall  be  divided  into  cells  centered 
at  xm,  >n,  Zp.  For  each  such  cell,  a  memory  and  processing  capability 
shall  be  provided.  Since  the  customary  echo-ranging  range/time  equation 
does  not  apply,  we  must  compute,  for  each  sensor /source/antenna -beam 
combination,  the  appropriate  range/time  relationship.  Figure  6  shows 
one  possible  technique.  For  ease  of  explanation,  a  horizontal  plane  is 
considered  in  this  figure.  Positions  Pj,  P^  and  ranges  Rj^,  Rja  +  RgJc 
are  known,  as  well  as  the  approximate  beam  angle  A.  It  is  desired  to 
classify  the  reported  R^a  +  Ralt  according  to  the  Xjj,  yn  cell(s)  it  may 
fit  into;  thus,  the  following  computations  are  needed: 


xm  "  Rak  008  A  +  *k 


y_  *  R  sin  A  +  y. 


Figure  6.  Range/time  relationship 


Clearly  must  be  found  as  a  function  of  Rak  +  R-a.  Constructing 

lines  B  and  D,  we  can  write  J 


B  sin  A  +  D  cos  A  =  x, 

k 


(33) 


B  cos  A 


D  sin  A  » 


(34) 


These  two  equations  in  B  and  D  are  readily  solved  for  each  beam  of  each 
sensor /source  combination*  Having  solved  for  B  and  D,  proceed  as 
follows: 

(a)  Divide  the  square  of  B  by  D  +  R^a  +  RaJt,  and  add  D  +  R^a  +  RaJt 
to  the  quotient.  Label  the  result  E: 


ja 


+  Rak  +  D 


+  R,  +  R 


^  • 

ja 


ak 


+  D 
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(b)  Show  that  E  =  2Rja  by  considering  that 


B2  +  (Rak  +  D)2  =  Rja  i 
82  =  Rja  -  (Rak  +  D)2  i 


82  "  (Rja  +  Rak  +  DHRja  '  Rak  "  D)  * 


- - - - -  =  R.  -  R  ,  -  D 

R.  +  R  .  +  D  na  ak 

3a  ak  J 


+  D  »  2R. 

3® 


Thus,  E  =  2Rja,  so  that  one  half  of  this  result  can  be  subtracted  from 
Rja  +  Rak  to  obtain  the  desired  value  of  Rak* 

When  the  reports  are  sorted  into  their  appropriate  xm,  yR,  z  cells, 
some  sorting  still  remains  to  be  done;  some  reports  will  certainly  fit 
into  more  than  one  cell,  because  of  sensor  beam  overlap,  coarseness  of 
the  sensor  beams,  and  uncertainty  of  beam  orientation.  In  addition,  it 
is  antipicated  that  more  than  one  report  in  a  given  source/sensor/beam 
combination  may  fit  into  a  given  cell,  this  will  certainly  happen  if  we 
wish  to  exploit  the  major  advantage  of  the  system— the  reliance  on  time 
to  determine  target  coordinates.  This  advantage  allows  us  to  achieve  a 
system  accuracy  that  greatly  exceeds  the  accuracy  of  its  components, 
thus  minimizing  the  number  of  range  cells  required  in  each  sensor. 


8.  DOPPLER  FREQUENCY  SHIFT  DUE  TO  MOTION  OF  SOURCE 

If  the  source  is  moving  with  velocity  vector  V,  the  envelope  of 
Doppler  effects  comprises  two  spheres  with  their  diameters  lying  on  the 
velocity  vector  and  their  common  point  of  tangency  at  the  source  (see 
fig.  7).  Each  sphere  has  diameter  D  =  V  and  represents  positive  or 
negative  Doppler,  respectively,  as  the  velocity  vector  points  toward  or 
away  from  it.  Diameter  D  is  assumed  small  compared  to  the  pulse 
repetition  frequency  (prf)  of  the  source.  Each  sensor,  to  achieve 
coherent  reception  of  echoes,  will  have  the  phase  of  its  local 
oscillator  adjusted,  pulse  by  pulse,  to  match  the  perceived  phase  of  the 


source  illumination  as  seen  via  the  direct  source-to-sensor  path.  The 
rate  at  which  this  phase  adjustment  occurs  represents  the  component  of 
source  Doppler  seen  by  the  sensor,  and  the  uniform  phase  progression 
experienced  by  the  local  oscillator  is  equivalent  to  a  frequency  shift 
that  makes  the  effective  frequency  of  the  local  oscillator  equal  to  the 
Doppler-shifted  frequency  of  the  source  as  seen  by  the  sensor. 


Figure  7.  Doppler  shifts  due  to  motion  of 
source . 


The  effect  of  this  circumstance  on  the  source  Doppler  picture  can  be 
visualized  as  follows: 

(a)  Find  the  "puncture  point"  where  the  source -to-sensor  line  of 
sight  pierces  the  source  Doppler  envelope. 

(b)  Construct  a  hemisphere,  centered  at  the  source,  with  its  base 
on  the  zero  Doppler  plane,  through  this  puncture  point. 


(c)  To  find  the  source  Doppler  as  seen  by  the  sensor  at  any  point, 
project  a  ray  from  the  source  to  that  point,  and  evaluate  the  vector 
length  lying  between  the  hemisphere  and  the  original  Doppler  envelope, 
on  this  ray. 

Doppler  frequency  shift  due  to  motion  of  the  source  can  be 
completely  characterized  at  the  system  level  as  follows: 

For  each  sensor,  k,  compute  the  direction  cosine  of  Rjk  as 


008  xjk  =  -rT 


yk  -  yj 

cos  =  “V" 


Z,  -  Z  . 

k  i 


008  zjk  “  ~ rT 


The  (unvnown)  components  of  the  source  Doppler  vector  are  D^x, 
Djy'  Djz* 

-The  Doppler  frequency  Fjk  observed  at  sensor  k  due  to  the  motion  of 
sou.wa  j  may  be  expressed  as 


Fjk  -  Djx  008  *jk  +  Djy  008  yjk  +  Djz  008  zj 


Putting  equations  (35)  in  the  matrix  format,  as  follows. 


cos  Xjl  cos  yjX  cos  Zjl  Djx 


cos  xjk  cos  yjk  cos  zjk  Djyi 


COS  Xj 


cos  y.,  cos  z 


jn  uj« 


•*  '4  *  *#■  *  . 4 


reveals  the  familiar  overcons trained  equations ,  and  we  proceed 
accordingly. 


The  implied  constraint  on  the  system  here  is  that  the  pulse  rate  of 
the  source  must  be  greater  than  twice  the  maximum  Doppler  because  of 
source  motion.  The  locus  of  rays  that  yields  zero  Doppler  defines  a 
cone  with  the  vertex  as  the  source  and  the  axis  along  the  source 
velocity  vector;  the  sensor  lies  on  the  surface  of  the  cone.  Similarly, 
any  locus  of  rays  that  yields  a  constant,  nonzero  Doppler  defines  a  cone 
with  the  same  vertex  and  axis,  differing  only  in  its  apex  angle. 

Thus,  with  a  high-altitude  source  and  a  sensor  located  on  a  planar 
terrain,  elliptical  (conic  section,  in  general)  isodops  with  the  zero 
isodop  passing  through  the  sensor  can  be  constructed;  however,  the 
sensor  will  usually  lie  on  a  major  axis  of  its  isodop.  The  major  axes 
of  all  isodops  for  all  source/sensor  combinations  will  lie  on  the 
vertical  projection  of  the  source  velocity  vector. 

All  the  foregoing  leads  to  the  conclusion  that  some  sensor 
processors  will  not  be  operating  in  the  best  possible  circumstances; 
this  is  especially  so  because  the  system  design  has  been  constrained  so 
that  no  sensor  has  any  of  the  analysis  developed  here  available  to  it. 

A  numerical  measure  of  the  source  Doppler  as  seen  by  each  sensor  can 
be  readily  generated,  and  these  numbers  can  be  transmitted  to  the 
central  processor  via  the  data  links  from  each  sensor  to  the  central 
processor.  The  carrier  frequency  of  these  data  links  can  serve  as  the 
frequency  reference  for  the  source  Doppler  measurements. 

With  the  source  Doppler  and  the  source -to-sensor  Doppler  available 
in  every  xn,  yu,  zp  cell,  their  difference  from  the  Doppler  as  reported 
on  each  alarm  from  every  sensor  can  be  constructed,  thus  eliminating 
Doppler  effects  due  to  source  motion.  After  the  alarm  reports  in  any 
given  cell  are  sorted,  the  velocity  vector  for  each  target  in  that  cell 
can  be  computed,  using  the  direction  cosines  of  the  cell-to-sensor  rays 
and  computing  as  was  done  in  determining  the  source  velocity  vector. 


9.  DISCUSSION 

It  is  evident  that  the  measurable  quantities  in  a  distributed 
surveillance  system  may  be  regarded  as  projections  upon  a  set  of  vectors 
that  is  established  by  the  relative  positions  of  the  sources  and/or 
sensors.  Hie  dimensionality  of  this  vector  space  grows  as  we  relinquish 
control  (or  absolute  knowledge)  of  the  relative  timing  of  events  within 
the  system.  It  is  interesting  to  note  that  lack  of  an  absolute  time 
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reference  puts  us  into  the  spacetime  geometry  required  by  special 
relativity's  discard  of  the  concept  of  simultaneity. 

Considering  the  spacelike  coordinates  alone,  the  relative 
independence  among  the  measurements  made  on  any  target  will  be 
proportional  to  the  size  of  the  solid  angle  subtended  by  the  source/ 
sensor  field  as  viewed  from  the  target.  Systems  having  a  small  solid 
angle  are  called  "skew  angular,"  since  all  their  basis  vectors  are 
skewed  from  the  desired  mutually  orthogonal  directions  into  a  single 
small  angular  spread. 

The  horrors  of  numerical  operations  in  highly  skew  angular 
coordinates  are  well-established3  and  can  serve  as  a  guide  to 
deployment,  with  the  proviso  that  the  total  spacetime  interval  (rather 
than  the  spacelike  coordinates  alone)  will  determine  the  degree  of 
difficulty  encountered  in  computation.  ttie  chief  computation  load, 
aside  from  each  sensor's  Doppler  filtering,  will  reside  in  the  target¬ 
sorting  algorithms.  Purely  incoherent  processing  of  sensor  reports  was 
the  only  technique  required  to  establish  complete  analysis  of  system 
geometry;  however,  it  should  be  clear  that  each  sensor  could  report  both 
amplitude  and  phase  of  target  Doppler,  allowing  further  coherent 
processing  at  the  central  computer. 

An  aspect  of  distributed  systems  whose  importance  cannot  be 
overemphasized  is  their  "distributed"  redundance  and  the  resultant  gain 
in  reliability  (or  operability)  that  they  provide.  The  Gaussian  least- 
squares  algorithm,  if  pursued  to  an  assessment  of  the  error  vector,  may 
be  regarded  as  an  introspective  process  that  has  miniscule  probability 
of  suffering  two  compensating  errors  that  will  produce  wrong  answers 
accompanied  by  an  acceptable  error  vector.  And  it  matters  not  where  in 
the  total  system  the  errors  occur.  In  the  detailed  discussion  of  the 
target-sorting  problem,  the  possibility  of  using  the  compatibility 
condition  was  not  mentioned,  for  lack  of  time  and  space.  Briefly,  we 
may  state  that  a  given  combination  of  target  reports  must  be  orthogonal 
to  all  solutions  of  the  adjoint  homogeneous  equations.  An  algorithm 
based  on  this  condition  may  provide  a  rapid  sorting  technique. 


10.  CONCLUSION 

The  geometry  of  distributed  surveillance  systems  is  amenable  to 
straightforward  calculation  in  linear  vector  spaces  of  no  more  than  five 
dimensions.  The  dimensionality  of  the  vector  space  depends  inversely  on 
the  knowledge  (or  control)  of  relative  timing  among  the  elements  of  the 


3 Cornelius  Lanczos,  Applied  Analysis,  Prentice-Hall  (1957) 


system.  Solutions  are  available  both  for  position  and  velocity,  and 
angular  resolution  at  sensors  is  required  only  to  ease  the  target¬ 
sorting  computational  load.  The  system  will  fail  soft,  and  immediately 
make  the  user  aware  of  all  but  a  most  improbable  set  of  errors  when  it 
fails.  The  systems  allow  a  remarkable  degree  of  autonomy  to  all 
elements,  and  retain  the  possibility  of  fully  coherent  processing  of  all 
sensor  signals.  Considerable  further  development  of  the  concepts 
addressed  herein  will  be  required  to  establish  the  total  feasibility  of 
any  distributed  surveillance  system,  but  the  possible  payoffs  are 

(a)  totally  unattended,  solar-powered  sensors,  emplaced  by 
helicopter  on  otherwise  unattainable  high  ground,  completely  passive 
except  for  a  low-probability-of-intercept  data  link,  or 

(b)  utterly  passive  tanks,  all  sharing  the  same  set  of  expendable 
sources,  with  each  tank  having  complete  knowledge  of  the  tactical 
situation  in  coordinates  centered  at  its  own  location. 
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